A Transcriptomic Classifier Model Identifies High-Risk Endotypes in a Prospective Study of Sepsis in Uganda.

OBJECTIVES: In high-income countries (HICs), sepsis endotypes defined by distinct pathobiological mechanisms, mortality risks, and responses to corticosteroid treatment have been identified using blood transcriptomics. The generalizability of these endotypes to low-income and middle-income countries (LMICs), where the global sepsis burden is concentrated, is unknown. We sought to determine the prevalence, prognostic relevance, and immunopathological features of HIC-derived transcriptomic sepsis endotypes in sub-Saharan Africa. DESIGN: Prospective cohort study. SETTING: Public referral hospital in Uganda. PATIENTS: Adults ( n = 128) hospitalized with suspected sepsis. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: Using whole-blood RNA sequencing data, we applied 19-gene and 7-gene classifiers derived and validated in HICs (SepstratifieR) to assign patients to one of three sepsis response signatures (SRS). The 19-gene classifier assigned 30 (23.4%), 92 (71.9%), and 6 (4.7%) patients to SRS-1, SRS-2, and SRS-3, respectively, the latter of which is designed to capture individuals transcriptionally closest to health. SRS-1 was defined biologically by proinflammatory innate immune activation and suppressed natural killer-cell, T-cell, and B-cell immunity, whereas SRS-2 was characterized by dampened innate immune activation, preserved lymphocyte immunity, and suppressed transcriptional responses to corticosteroids. Patients assigned to SRS-1 were predominantly (80.0% [24/30]) persons living with HIV with advanced immunosuppression and frequent tuberculosis. Mortality at 30-days differed significantly by endotype and was highest (48.1%) in SRS-1. Agreement between 19-gene and 7-gene SRS assignments was poor (Cohen's kappa 0.11). Patient stratification was suboptimal using the 7-gene classifier with 15.1% (8/53) of individuals assigned to SRS-3 deceased at 30-days. CONCLUSIONS: Sepsis endotypes derived in HICs share biological and clinical features with those identified in sub-Saharan Africa, with major differences in host-pathogen profiles. Our findings highlight the importance of context-specific sepsis endotyping, the generalizability of conserved biological signatures of critical illness across disparate settings, and opportunities to develop more pathobiologically informed sepsis treatment strategies in LMICs.

Development of a Novel Clinicomolecular Risk Index to Enhance Mortality Prediction and Immunological Stratification of Adults Hospitalized with Sepsis in Sub-Saharan Africa: A Pilot Study from Uganda.

The global burden of sepsis is concentrated in sub-Saharan Africa (SSA), where epidemic HIV and unique pathogen diversity challenge the effective management of severe infections. In this context, patient stratification based on biomarkers of a dysregulated host response may identify subgroups more likely to respond to targeted immunomodulatory therapeutics. In a prospective cohort of adults hospitalized with suspected sepsis in Uganda, we applied machine learning methods to develop a prediction model for 30-day mortality that integrates physiology-based risk scores with soluble biomarkers reflective of key domains of sepsis immunopathology. After model evaluation and internal validation, whole-blood RNA sequencing data were analyzed to compare biological pathway enrichment and inferred immune cell profiles between patients assigned differential model-based risks of mortality. Of 260 eligible adults (median age, 32 years; interquartile range, 26-43 years; 59.2% female, 53.9% living with HIV), 62 (23.8%) died by 30 days after hospital discharge. Among 14 biomarkers, soluble tumor necrosis factor receptor 1 (sTNFR1) and angiopoietin 2 (Ang-2) demonstrated the greatest importance for mortality prediction in machine learning models. A clinicomolecular model integrating sTNFR1 and Ang-2 with the Universal Vital Assessment (UVA) risk score optimized 30-day mortality prediction across multiple performance metrics. Patients assigned to the high-risk, UVA-based clinicomolecular subgroup exhibited a transcriptional profile defined by proinflammatory innate immune and necroptotic pathway activation, T-cell exhaustion, and expansion of key immune cell subsets including regulatory and gamma-delta T cells. Clinicomolecular stratification of adults with suspected sepsis in Uganda enhanced 30-day mortality prediction and identified a high-risk subgroup with a therapeutically targetable immunological profile. Further studies are needed to advance pathobiologically informed sepsis management in SSA.

Brief Report: Detection of Urine Lipoarabinomannan Is Associated With Proinflammatory Innate Immune Activation, Impaired Host Defense, and Organ Dysfunction in Adults With Severe HIV-Associated Tuberculosis in Uganda.

BACKGROUND: The immunopathology of disseminated HIV-associated tuberculosis (HIV/TB), a leading cause of critical illness and death among persons living with HIV in sub-Saharan Africa, is incompletely understood. Reflective of hematogenously disseminated TB, detection of lipoarabinomannan (LAM) in urine is associated with greater bacillary burden and poor outcomes in adults with HIV/TB. METHODS: We determined the relationship between detection of urine TB-LAM, organ dysfunction, and host immune responses in a prospective cohort of adults hospitalized with severe HIV/TB in Uganda. Generalized additive models were used to analyze the association between urine TB-LAM grade and concentrations of 14 soluble immune mediators. Whole-blood RNA-sequencing data were used to compare transcriptional profiles between patients with high- vs. low-grade TB-LAM results. RESULTS: Among 157 hospitalized persons living with HIV, 40 (25.5%) had positive urine TB-LAM testing. Higher TB-LAM grade was associated with more severe physiologic derangement, organ dysfunction, and shock. Adjusted generalized additive models showed that higher TB-LAM grade was significantly associated with higher concentrations of mediators reflecting proinflammatory innate and T-cell activation and chemotaxis (IL-8, MIF, MIP-1ß/CCL4, and sIL-2Ra/sCD25). Transcriptionally, patients with higher TB-LAM grades demonstrated multifaceted impairment of antibacterial defense including reduced expression of genes encoding cytotoxic and autophagy-related proteins and impaired cross-talk between innate and cell-mediated immune effectors. CONCLUSIONS: Our findings add to emerging data suggesting pathobiological relationships between LAM, TB dissemination, innate cell activation, and evasion of host immunity in severe HIV/TB. Further translational studies are needed to elucidate the role for immunomodulatory therapies, in addition to optimized anti-TB treatment, in this often critically ill population.

HIV infection drives pro-inflammatory immunothrombotic pathway activation and organ dysfunction among adults with sepsis in Uganda.

BACKGROUND: The global burden of sepsis is concentrated in high HIV-burden settings in sub-Saharan Africa (SSA). Despite this, little is known about the immunopathology of sepsis in persons with HIV (PWH) in the region. We sought to determine the influence of HIV on host immune responses and organ dysfunction among adults hospitalized with suspected sepsis in Uganda. DESIGN: Prospective cohort study. METHODS: We compared organ dysfunction and 30-day outcome profiles of PWH and those without HIV. We quantified 14 soluble immune mediators, reflective of key domains of sepsis immunopathology, and performed whole-blood RNA-sequencing on samples from a subset of patients. We used propensity score methods to match PWH and those without HIV by demographics, illness duration, and clinical severity, and compared immune mediator concentrations and gene expression profiles across propensity score-matched groups. RESULTS: Among 299 patients, 157 (52.5%) were PWH (clinical stage 3 or 4 in 80.3%, 67.7% with known HIV on antiretroviral therapy). PWH presented with more severe physiologic derangement and shock, and had higher 30-day mortality (34.5% vs. 10.2%; P  < 0.001). Across propensity score-matched groups, PWH exhibited greater pro-inflammatory immune activation, including upregulation of interleukin (IL)-6, IL-8, IL-15, IL-17 and HMGB1 signaling, with concomitant T-cell exhaustion, prothrombotic pathway activation, and angiopoeitin-2-related endothelial dysfunction. CONCLUSIONS: Sepsis-related organ dysfunction and mortality in Uganda disproportionately affect PWH, who demonstrate exaggerated activation of multiple immunothrombotic and metabolic pathways implicated in sepsis pathogenesis. Further investigations are needed to refine understanding of sepsis immunopathology in PWH, particularly mechanisms amenable to therapeutic manipulation.

Dynamics of SARS-CoV-2 variants characterized during different COVID-19 waves in Mali.

Background: The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) variants may have contributed to prolonging the pandemic, and increasing morbidity and mortality related to coronavirus disease 2019 (COVID-19). This article describes the dynamics of circulating SARS-CoV-2 variants identified during the different COVID-19 waves in Mali between April and October 2021. Methods: The respiratory SARS-CoV-2 complete spike (S) gene from positive samples was sequenced. Generated sequences were aligned by Variant Reporter v3.0 using the Wuhan-1 strain as the reference. Mutations were noted using the GISAID and Nextclade platforms. Results: Of 16,797 nasopharyngeal swab samples tested, 6.0% (1008/16,797) tested positive for SARS-CoV-2 on quantitative reverse transcription polymerase chain reaction. Of these, 16.07% (162/1008) had a cycle threshold value ≤28 and were amplified and sequenced. The complete S gene sequence was recovered from 80 of 162 (49.8%) samples. Seven distinct variants were identified: Delta (62.5%), Alpha (1.2%), Beta (1.2%), Eta (30.0%), 20B (2.5%), 19B (1.2%) and 20A (1.2%). Conclusions and perspectives: Several SARS-CoV-2 variants were present during the COVID-19 waves in Mali between April and October 2021. The continued emergence of new variants highlights the need to strengthen local real-time sequencing capacity and genomic surveillance for better and coordinated national responses to SARS-CoV-2.

Detection of antibodies to Anaplasma phagocytophilum and Babesia microti using linear peptides.

Anaplasma phagocytophilum and Babesia microti are emerging tick-borne pathogens in the United States. Although active infection is typically diagnosed by direct diagnostic tests, such as blood smear or polymerase chain reaction assay, serologic assays can be helpful to identify past infections, and the use of acute plus convalescent testing can potentially identify recent infections. We employed a peptide array to select sets of linear peptides for serologic diagnosis of infections with A. phagocytophilum and B. microti. Three optimal peptides were selected for each agent based on their performance with clinical specimens. All three A. phagocytophilum peptides were located within the conserved fragments of the MSP2 antigen. Two B. microti peptides were located in the N terminus of the SA-1 antigen; the third was in the BMN 1-17 antigen. We found that these peptides can be a useful tool for detection of antibody reactivity to both of these pathogens.

Virome of Ixodes ricinus, Dermacentor reticulatus, and Haemaphysalis concinna Ticks from Croatia.

Tick-borne diseases are a serious threat to both public and veterinary health. In this study, we used high-throughput sequencing to characterize the virome of three tick species implicated in the spread of vector-borne disease throughout Croatia. Ten viruses were identified, including seven potential novel species within the viral families Flaviviridae, Nyamiviridae, Rhabdoviridae, Peribunyaviridae, Phenuiviridae, and Nairoviridae.

Multidimensional analysis of the host response reveals prognostic and pathogen-driven immune subtypes among adults with sepsis in Uganda.

BACKGROUND: The global burden of sepsis is concentrated in sub-Saharan Africa, where severe infections disproportionately affect young, HIV-infected adults and high-burden pathogens are unique. In this context, poor understanding of sepsis immunopathology represents a crucial barrier to development of locally-effective treatment strategies. We sought to determine inter-individual immunologic heterogeneity among adults hospitalized with sepsis in a sub-Saharan African setting, and characterize associations between immune subtypes, infecting pathogens, and clinical outcomes. METHODS: Among a prospective observational cohort of 288 adults hospitalized with suspected sepsis in Uganda, we applied machine learning methods to 14 soluble host immune mediators, reflective of key domains of sepsis immunopathology (innate and adaptive immune activation, endothelial dysfunction, fibrinolysis), to identify immune subtypes in randomly-split discovery (N = 201) and internal validation (N = 87) sub-cohorts. In parallel, we applied similar methods to whole-blood RNA-sequencing data from a consecutive subset of patients (N = 128) to identify transcriptional subtypes, which we characterized using biological pathway and immune cell-type deconvolution analyses. RESULTS: Unsupervised clustering consistently identified two immune subtypes defined by differential activation of pro-inflammatory innate and adaptive immune pathways, with transcriptional evidence of concomitant CD56(-)/CD16( +) NK-cell expansion, T-cell exhaustion, and oxidative-stress and hypoxia-induced metabolic and cell-cycle reprogramming in the hyperinflammatory subtype. Immune subtypes defined by greater pro-inflammatory immune activation, T-cell exhaustion, and metabolic reprogramming were consistently associated with a high-prevalence of severe and often disseminated HIV-associated tuberculosis, as well as more extensive organ dysfunction, worse functional outcomes, and higher 30-day mortality. CONCLUSIONS: Our results highlight unique host- and pathogen-driven features of sepsis immunopathology in sub-Saharan Africa, including the importance of severe HIV-associated tuberculosis, and reinforce the need to develop more biologically-informed treatment strategies in the region, particularly those incorporating immunomodulation.

Viral capture sequencing detects unexpected viruses in the cerebrospinal fluid of adults with meningitis.

OBJECTIVES: Many patients with meningitis have no aetiology identified leading to unnecessary antimicrobials and prolonged hospitalisation. We used viral capture sequencing to identify possible pathogenic viruses in adults with community-acquired meningitis. METHODS: Cerebrospinal fluid (CSF) from 73 patients was tested by VirCapSeq-VERT, a probe set designed to capture viral targets using high throughput sequencing. Patients were categorised as suspected viral meningitis - CSF pleocytosis, no pathogen identified (n = 38), proven viral meningitis - CSF pleocytosis with a pathogen identified (n = 15) or not meningitis - no CSF pleocytosis (n = 20). RESULTS: VirCapSeq-VERT detected virus in the CSF of 16/38 (42%) of those with suspected viral meningitis, including twelve individual viruses. A potentially clinically relevant virus was detected in 9/16 (56%). Unexpectedly Toscana virus, rotavirus and Saffold virus were detected and assessed to be potential causative agents. CONCLUSION: VirCapSeq-VERT increases the probability of detecting a virus. Using this agnostic approach we identified Toscana virus and, for the first time in adults, rotavirus and Saffold virus, as potential causative agents in adult meningitis. Further work is needed to determine the prevalence of atypical viral candidates as well as the clinical impact of using sequencing methods in real time. This knowledge can help to reduce antimicrobial use and hospitalisations leading to both patient and health system benefits.

SARS-CoV-2 Sequence Analysis during COVID-19 Case Surge, Liberia, 2021.

In June 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases surged in Liberia. SARS-CoV-2 sequences from patients hospitalized during March-July 2021 revealed the Delta variant was in Liberia in early March and was dominant in June, irrespective of geography. Mutations and deletions suggest multiple SARS-CoV-2 Delta variant introductions.

Development of a capture sequencing assay for enhanced detection and genotyping of tick-borne pathogens.

Inadequate sensitivity has been the primary limitation for implementing high-throughput sequencing for studies of tick-borne agents. Here we describe the development of TBDCapSeq, a sequencing assay that uses hybridization capture probes that cover the complete genomes of the eleven most common tick-borne agents found in the United States. The probes are used for solution-based capture and enrichment of pathogen nucleic acid followed by high-throughput sequencing. We evaluated the performance of TBDCapSeq to surveil samples that included human whole blood, mouse tissues, and field-collected ticks. For Borrelia burgdorferi and Babesia microti, the sensitivity of TBDCapSeq was comparable and occasionally exceeded the performance of agent-specific quantitative PCR and resulted in 25 to > 10,000-fold increase in pathogen reads when compared to standard unbiased sequencing. TBDCapSeq also enabled genome analyses directly within vertebrate and tick hosts. The implementation of TBDCapSeq could have major impact in studies of tick-borne pathogens by improving detection and facilitating genomic research that was previously unachievable with standard sequencing approaches.

Novel quaranjavirus and other viral sequences identified from ticks parasitizing hunted wildlife in Trinidad and Tobago.

Hunters are at a higher risk for exposure to zoonotic pathogens due to their close interactions with wildlife and arthropod vectors. In this study, high throughput sequencing was used to explore the viromes of two tick species, Amblyomma dissimile and Haemaphysalis juxtakochi, removed from hunted wildlife in Trinidad and Tobago. We identified sequences from 3 new viral species, from the viral families Orthomyxoviridae, Chuviridae and Tetraviridae in A. dissimile.

A metagenomic examination of the pathobiome of the invasive tick species, Haemaphysalis longicornis, collected from a New York City borough, USA.

Haemaphysalis longicornis, the Asian longhorned tick, is an invasive tick species that has spread rapidly across the northeastern and southeastern regions of the United States in recent years. This invasive pest species, known to transmit several tick-borne pathogens in its native range, is a potential threat to wildlife, livestock, domestic animals, and humans. Questing larval (n = 25), nymph (n = 10), and adult (n = 123), along with host-derived adult (n = 25) H. longicornis ticks were collected from various locations on Staten Island, NY. The pathobiome of each specimen was examined using two different high throughput sequencing approaches, virus enrichment and shotgun metagenomics. An average of 45,828,061 total reads per sample were recovered from the virus enriched samples and an average of 11,381,144 total reads per sample were obtained using shotgun metagenomics. Aside from endogenous viral sequences, no viruses were identified through either approach. Through shotgun metagenomics, Coxiella-like bacteria, Legionella, Sphingomonas, and other bacterial species were recovered. The Coxiella-like agent was ubiquitous and present at high abundances in all samples, suggesting it may be an endosymbiont. The other bacterial agents are not known to be transmitted by ticks. From these analyses, H. longicornis do not appear to host any endemic human tick-borne pathogens in the New York City region.

Impact of CMV Reactivation, Treatment Approaches, and Immune Reconstitution in a Nonmyeloablative Tolerance Induction Protocol in Cynomolgus Macaques.

BACKGROUND: Cytomegalovirus (CMV) infection is a serious complication in immunosuppressed patients, specifically transplant recipients. Here, we describe the development and use of an assay to monitor the incidence and treatment of CMV viremia in a Cynomolgus macaque model of bone marrow transplantation (BMT) for tolerance induction. We address the correlation between the course of viremia and immune reconstitution. METHODS: Twenty-one animals received a nonmyeloablative conditioning regimen. Seven received cyclosporine A for 28 days and 14 received rapamycin. A CMV polymerase chain reaction assay was developed and run twice per week to monitor viremia. Nineteen recipients were CMV seropositive before BMT. Immune reconstitution was monitored through flow cytometry and CMV viremia was tracked via quantitative polymerase chain reaction. RESULTS: Recipients developed CMV viremia during the first month post-BMT. Two animals developed uncontrollable CMV disease. CMV reactivation occurred earlier in cyclosporine A-treated animals compared with those receiving rapamycin. Post-BMT, T-cell counts remained significantly lower compared with pretransplant levels until CMV reactivation, at which point they increased during the viremic phase and approached pretransplant levels 3 months post-BMT. Management of CMV required treatment before viremia reached 10 000 copies/mL; otherwise clinical symptoms were observed. High doses of ganciclovir resolved the viremia, which could subsequently be controlled with valganciclovir. CONCLUSIONS: We developed an assay to monitor CMV in Cynomolgus macaques. CMV reactivation occurred in 100% of seropositive animals in this model. Rapamycin delayed CMV reactivation and ganciclovir treatment was effective at high doses. As in humans, CD8 T cells proliferated during CMV viremia.

Identification of immunoreactive linear epitopes of Borrelia miyamotoi.

Borrelia miyamotoi is an emerging tick-borne spirochete transmitted by ixodid ticks. Current serologic assays for B. miyamotoi are impacted by genetic similarities to other Borrelia and limited understanding of optimal antigenic targets. In this study, we employed the TBD-Serochip, a peptide array platform, to identify new linear targets for serologic detection of B. miyamotoi. We examined a wide range of suspected B. miyamotoi antigens and identified 352 IgM and 91 IgG reactive peptides, with the majority mapping to variable membrane proteins. These included peptides within conserved fragments of variable membrane proteins that may have greater potential for differential diagnosis. We also identified reactive regions on FlaB, and demonstrate crossreactivity of B. burgdorferi s.l. C6 with a B. miyamotoi C6-like peptide. The panel of linear peptides identified in this study can be used to enhance serodiagnosis of B. miyamotoi.

Viral Diversity of Tick Species Parasitizing Cattle and Dogs in Trinidad and Tobago.

Ticks are vectors of a wide variety of pathogens that are implicated in mild to severe disease in humans and other animals. Nonetheless, the full range of tick-borne pathogens is unknown. Viruses, in particular, have been neglected in discovery efforts targeting tick-borne agents. High throughput sequencing was used to characterize the virome of 638 ticks, including Rhipicephalus microplus (n = 320), Rhipicephalus sanguineus (n = 300), and Amblyomma ovale (n = 18) collected throughout Trinidad and Tobago in 2017 and 2018. Sequences representing nine viruses were identified, including five novel species within Tymovirales, Bunyavirales, Chuviridae, Rhabdoviridae, and Flaviviridae. Thereafter the frequency of detection of viral sequences in individual tick species was investigated.

Virome capture sequencing does not identify active viral infection in unicentric and idiopathic multicentric Castleman disease.

Castleman disease (CD) describes a spectrum of heterogeneous disorders defined by characteristic lymph node histopathology. Enlarged lymph nodes demonstrating CD histopathology can occur in isolation (unicentric CD; UCD) sometimes accompanied by mild symptoms, or at multiple sites (multicentric CD, MCD) with systemic inflammation and cytokine-driven multi-organ dysfunction. The discovery that Kaposi sarcoma herpesvirus/human herpesvirus (HHV)-8 drives MCD in a subset of patients has led to the hypotheses that UCD and MCD patients with negative HHV-8 testing by conventional methods may represent false negatives, or that these cases are driven by another virus, known or unknown. To investigate these hypotheses, the virome capture sequencing for vertebrate viruses (VirCapSeq-VERT) platform was employed to detect RNA transcripts from known and novel viruses in fresh frozen lymph node tissue from CD patients (12 UCD, 11 HHV-8-negative MCD [idiopathic MCD; iMCD], and two HHV-8-positive MCD) and related diseases (three T cell lymphoma and three Hodgkin lymphoma). This assay detected HHV-8 in both HHV-8-positive cases; however, HHV-8 was not found in clinically HHV-8-negative iMCD or UCD cases. Additionally, no novel viruses were discovered, and no single known virus was detected with apparent association to HHV-8-negative CD cases. Herpesviridae family members, notably including Epstein-Barr virus (EBV), were detected in 7 out of 12 UCD and 5 of 11 iMCD cases with apparent correlations with markers of disease severity in iMCD. Analysis of a separate cohort of archival formalin-fixed, paraffin-embedded lymph node tissue by In situ hybridization revealed significantly fewer EBV-positive cells in UCD and iMCD compared to tissue from HHV-8-positive MCD and EBV-associated lymphoproliferative disorder. In an additional cohort, quantitative testing for EBV by PCR in peripheral blood during disease flare did not detect systemic EBV viremia, suggesting detection lymph node tissue is due to occult, local reactivation in UCD and iMCD. This study confirms that HHV-8 is not present in UCD and iMCD patients. Further, it fails to establish a clear association between any single virus, novel or known, and CD in HHV-8-negative cases. Given that distinct forms of CD exist with viral and non-viral etiological drivers, CD should be considered a group of distinct and separate diseases with heterogeneous causes worthy of further study.

Characterization of Stool Virome in Children Newly Diagnosed With Moderate to Severe Ulcerative Colitis.

BACKGROUND: Viral infections have been suggested as possible triggers for the onset of ulcerative colitis (UC). METHODS: We employed VirCapSeq-Vert, a high-throughput sequencing virus capture platform, to examine the stool virome of children with newly diagnosed moderate to severe UC. We surveyed fecal samples collected at presentation, after symptom remission, and from a control group diagnosed with irritable bowel syndrome. RESULTS: Seventy subjects with UC (mean age 13 years, 45 had moderate symptoms, 25 had severe, 69 of 70 had a Mayo endoscopy subscore 2/3) were studied. We detected a wide range of animal viruses that were taxonomically classified into 12 viral families. A virus was present in 50% of fecal samples collected at presentation, 41% of samples collected after remission, and 40% of samples in our control group. The most frequently identified viruses were diet-based gyroviruses. The UC cohort had a significantly higher prevalence of anelloviruses compared with the control cohort. However, we did not identify a single virus that can be implicated in the onset of UC and did not find an association between UC disease severity and viral presence. CONCLUSION: Presence of virus in stool was not associated with the onset of pediatric UC.

Microbiome analysis of Ixodes scapularis ticks from New York and Connecticut.

We employed high throughput sequencing to survey the microbiomes of Ixodes scapularis collected in New York and Connecticut. We examined 197 individual I. scapularis adults and pools from 132 adults and 197 nymphs. We detected Borrelia burgdorferi sensu stricto in 56.3% of individual ticks, Anaplasma phagocytophilum in 10.6%, Borrelia miyamotoi in 5%, Babesia microti in 7.6%, and Powassan virus in 3.6%. We did not detect Borrelia mayonii, Ehrlichia muris eauclairensis, Bartonella spp. or pathogenic Babesia species other than B. microti. The most abundant bacterium (65%), and only rickettsial species identified, was the endosymbiont Rickettsia buchneri. A filarial nematode was found in 13.7% of adult ticks. Fourteen viruses were detected including South Bay virus (22%) and blacklegged tick phlebovirus 1 and 2 (73%). This study provides insight into the microbial diversity of I. scapularis in New York State and Connecticut.

Identification of Novel Viruses in Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis Ticks.

Ticks carry a wide range of known human and animal pathogens and are postulated to carry others with the potential to cause disease. Here we report a discovery effort wherein unbiased high-throughput sequencing was used to characterize the virome of 2,021 ticks, including Ixodes scapularis (n = 1,138), Amblyomma americanum (n = 720), and Dermacentor variabilis (n = 163), collected in New York, Connecticut, and Virginia in 2015 and 2016. We identified 33 viruses, including 24 putative novel viral species. The most frequently detected viruses were phylogenetically related to members of the Bunyaviridae and Rhabdoviridae families, as well as the recently proposed Chuviridae. Our work expands our understanding of tick viromes and underscores the high viral diversity that is present in ticks. IMPORTANCE The incidence of tick-borne disease is increasing, driven by rapid geographical expansion of ticks and the discovery of new tick-associated pathogens. The examination of the tick microbiome is essential in order to understand the relationship between microbes and their tick hosts and to facilitate the identification of new tick-borne pathogens. Genomic analyses using unbiased high-throughput sequencing platforms have proven valuable for investigations of tick bacterial diversity, but the examination of tick viromes has historically not been well explored. By performing a comprehensive virome analysis of the three primary tick species associated with human disease in the United States, we gained substantial insight into tick virome diversity and can begin to assess a potential role of these viruses in the tick life cycle.